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Abstract 

The jets observed to emanate from many compact accreting objects 
may arise from the twisting of a magnetic field threading a differentiaUy 
rotating accretion disk which acts to magneticaUy extract angular momen- 
tum and energy from the disk. Two main regimes have been discussed, 
hydromagnetic jets, which have a significant mass flux and have energy 
and angular momentum carried by both matter and electromagnetic held 
and, Poynting jets, where the mass flux is small and energy and angu- 
lar momentum are carried predominantly by the electromagnetic field. 
Here, we describe recent theoretical work on the formation of relativis- 
tic Poynting jets from magnetized accretion disks. Further, we describe 
new relativistic, fully-electromagnetic, particle-in-cell simulations of the 
formation of jets from accretion disks. Analog Z-pinch experiments may 
help to understand the origin of astrophysical jets. 

1. INTRODUCTION 

Powerful, highly-coUimated, oppositely directed jets are observed in active 
galaxies and quasars (see for example Bridle & Eilek 1984), and in old compact 
stars in binaries - the "microquasars" (Mirabel & Rodriguez 1994; Eikenberry et 
al. 1998). Further, highly coUimated emission line jets are seen in young stellar 
objects (Biihrke, Mundt, & Ray 1988). Different models have been put for- 
ward to explain astrophysical jets (Bisnovatyi-Kogan & Lovelace 2001). Recent 
observational and theoretical work favors models where twisting of an ordered 
magnetic field threading an accretion disk acts to magnetically accelerate the 
jets. Here, we discuss the origin of the relativistic jets observed in active galax- 
ies and quasars and in microquasars. We first discuss a theoretical model (§1), 
and then new results from relativistic particle-in-cell (PIC) simulations (§2). 

2. POYNTING JETS 

The powerful jets observed from active galaxies and quasars are probably 
not hydromagnetic outflows but rather Poynting flux dominated jets. The mo- 
tions of these jets measured by very long baseline interferometry correspond to 
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bulk Lorentz factors of F = 0(10) which is much larger than the Lorentz factor 
of the Keplerian disk velocity predicted for hydromagnetic outflows. Further- 
more, the low Faraday rotation measures observed for these jets at distances < 
kpc from the central object implies a very low plasma densities. Similar argu- 
ments indicate that the jets of microquasars are not hydromagnetic outflows but 
rather Poynting jets. Poynting jets have also been proposed to be the driving 
mechanism for gamma ray burst sources (Katz 1997). Theoretical studies have 
developed models for Poynting jets from accretion disks (Lovelace, Wang, & 
Sulkanen 1987; Lyndcn-Bell 1996; Romanova & Lovelace 1997; Lcvinson 1998; 
Li et al. 2001; Lovelace et al. 2002; and Lovelace & Romanova 2003). Station- 
ary non-relativistic Poynting flux dominated outflows were found by Romanova 
et al. (1998) and Ustyugova et al. (2000) in axisymmetric MHD simulations of 
the opening of magnetic loops threading a Keplerian disk. Here, we summma- 
rize a model for the formation of relativistic Poynting jets from a disk (Lovelace 
& Romanova 2003). 

Consider a dipole-like coronal magnetic field - such as that shown in the 
lower part of Figure la - threading a differentially rotating Keplerian accretion 
disk. The disk is perfectly conducting, high-density, and has a small accretion 
speed (^ Keplerian speed). The field may be generated in the disk by a dynamo 
action and released. Outside of the disk there is assumed to be a "coronal" or 
"force-free" plasma (pgE-f-J x B/c « 0, Gold & Hoyle 1960). We use cylindrical 
(r, 0, z) coordinates and consider axisymmetric field configurations. Thus the 
magnetic field has the form B = Bp -|- with Bp = B^t + B^z. We have 
Br = -{l/r)d'^/dz and B^ = {l/r)d^/dr. where ^{r,z) = rA^{r,z) is the 
flux function. 

Most of the azimuthal twist of a fleld line of the Poynting jet occurs 

along the jet from z = to Z{t?} as sketched in Figure la, where Z(t) is the 
axial location of the "head" of the jet. Along most of the distance z = Oto Z, the 
radius of the jet is a constant and = ^'(r) for Z >> ro, where ro is the radius 
of the 0-point of the magnetic field in the disk. Note that the function ^(r) 
is different from \l/(r, 0) which is the flux function profile on the disk surface. 
Hence r^d(f)/dz = rB^{r,z)/Bz{r,z). We take for simplicity Vz = dZ/dt = 
const. We determine Vz subsequently. In this case H{^>) = [r'^n{'ii) /Vz]Bz 
can be written as a function of 'f and d^ /dr. With H known, the relativistic 
Grad-Shafranov equation, 



can be solved (Lovelace & Romanova 2003). 

The quantity not determined by equation (1) is the velocity 14, or Lorentz 
factor r = 1/(1 — Ic^yl"^. This is determined by taking into account the 
balance of axial forces at the head of the jet: the electromagnetic pressure within 
the jet is balanced against the dynamic pressure of the external medium which 
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is assumed uniform with density pext- This gives (F^ — 1)'^ = Bq / {SttTZ"^ pextc^) , 
or for r > 1, 

where TZ = ro/vg ^ 1 and Vg = GM/c^, and Bq the magnetic field strength at 
the center of the disk. A necessary condition for the vahdity of this equation 
is that the axial speed of the counter-propagating fast magnetosonic wave (in 
the lab frame) be larger than Vz so that the jet is effectively 'subsonic' This 
value of r is of the order of the Lorentz factors of the expansion of parsec- 
scale extragalactic radio jets observed with very-long-baseline-interferometry 
(see, e.g., Zensus et al. 1998). This interpretation assumes that the radiating 
electrons (and/or positrons) are accelerated to high Lorentz factors (7 ~ 10^) at 
the jet front and move with a bulk Lorentz factor V relative to the observer. The 
luminosity of the +z Poynting jet is Ej = c rdrErB^/2 = cB^Te^/^r^/S ~ 
2.1 X lO46(Bo/lO3G)2(7e/lO)3/2(M/lO9M0)2 erg/s, where M is the mass of the 
black hole. 

For long time-scales, the Poynting jet is time-dependent due to the angular 
momentum it extracts from the inner disk (r < ro) which in turn causes ro 
to decrease with time (Lovelace et al. 2002). This loss of angular momentum 
leads to a "global magnetic instability" and collapse of the inner disk (Lovelace 
et al. 1994, 1997, 2002) and a corresponding outburst of energy in the jets 
from the two sides of the disk. Such outbursts may explain the flares of active 
galactic nuclei blazar sources (Romanova & Lovelace 1997; Levinson 1998) and 
the one-time outbursts of gamma ray burst sources (Katz 1997). 

3. RELATIVISTIC ELECTROMAGNETIC PIC SIMULATIONS 

We performed relativistic, fully electromagnetic, particle-in-cell simulations 
of the formation of jets from an acc;rction disk initially threaded by a dipole- 
like magnetic field. This was done using the code XOOPIC developed by Ver- 
boncoeur, Langdon, and Gladd (1995). Earlier, Gisler, Lovelace, and Norman 
(1989) studied jet formation for a monopole type field using the relativistic PIC 
code ISIS. The geometry of the initial configuration is shown in Figure lb. The 
computational region is a cylindrical "can," r = — Rm and z = — Z^, with 
outflow boundary conditions on the outer boundaries, and the potential and 
particle emission specified on the disk surface r = — Rm, z = 0. Equal fluxes 
of electrons and positrons are emitted so that the net emission is effectively 
space-charge-limited. About 10^ particles were used in the simulations reported 
here. The behavior of the lower half-space {z < 0) is expected to be a mirror 
image of the upper half-space. 

Figure 2 shows the formation of a relativistic jet. The gray scale indicates 
the logarithm of the density of electrons or positrons with 20 levels between the 
lightest (10^^) and darkest (4 x 10^^/m^). The lines are poloidal magnetic field 
lines Bp. The total B— field is shown in Figure 3. The computational region 



3 



has {Rm,Zm) = (50,100) m, the initial B— field is dipole-like with iJz(0,0) = 
-Bo = 28.3 G and an 0-point at (r, z) = (10,0) m, and the electric potential 
at the center of the disk is $o = —10^ V relative to the outer region of the 
disk. Initially, the computational region was filled with a distribution of equal 
densities of electrons and positrons with n±(0, 0) = 3 x lO^^/m^. Electrons and 
positrons are emitted with equal currents /± = 3 x 10^ A from both the inner 
and the outer portions of the disk as indicated in Figure lb with an axial speed 
much less than c. For a Keplerian disk with ro ^ rg, the scalings are $o ~ 
BQirorg)^^"^ , I ^ cBoTq and the jet power is ^ cBgr^^^ry^ . The calculations 
were done on a 64 x 128 grid stretched in both the r and z directions so as to 
give much higher spatial resolution at small r and small z. These simulations 
show the formation of a quasi-stationary, coUimated current-carrying jet. The 
Poynting flux power of the jet is Ej ^ 7x 10^^ W and the particle kinetic energy 
power is » 4.7 x 10^° W. The charge density of the electron flow is partially 
neutralized by the positron flow. Simlations are planned with the positrons 
replaced by ions. Scaled Z-pinch experiments configured as shown in Figure lb 
can allow further study of astrophysical jet formation. 
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J. Wick, and T. Womack for valuable assistance and discussions. This work 
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Figure 1: (a) Sketch of the magnetic field configuration of a Poynting jet from 
Lovelace and Romanova (2003). The bottom part of (a) shows the initial dipolc- 
like magnetic field threading the disk which rotates at the angular rate f2(r). The 
top part shows the jet at some time later when the head of the jet is at a distance 
Z{t). At the head of the jet there is force balance between electromagnetic 
stress of the jet and the ram pressure of the ambient medium of density pext- 
(b) Sketch of the initial conditions for the relativistic PIC simulations of jet 
formation from an accretion disk. 

Mirabel, I.F., & Rodriguez, L.F. 1994 Nature, 371, 46 
Romanova M.M., & Lovelace R.V.E. 1997, ApJ, 475, 97 

Romanova, M.M., Ustyugova, G.V, Koldoba, A.V, Chechetkin, VM., & Lovelace, 
R.V.E. 1998, ApJ, 500, 703 

Ustyugova, G.V, Lovelace, R.V.E. , Romanova, M.M., Li, H., & Colgate, S.A. 
2000 ApJ, 541, L21 

Verboncoeur, J. P., Langdon, A.B., & Gladd, N.T. 1995, Comp. Phys. Comm., 
87, 199 

Zensus, J. A., Taylor, G.B., & Wrobcl, J.M. (cds.) 1998, Radio Emission from 
Galactic and Extragalactic Compact Sources, lAU Colloquium 164, (Ast. 
Soc. of the Pacific) 



5 



Figure 2: Relativistic PIC simulations of the formation of a jet from a rotating 
disk, (a) -(c) give snapshots at times (1, 2, 3) x 10~^ s, and (d) is at f = 10~^ s. 
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Figure 3: Three dimensional magnetic field lines originating from the disk at 
r = 1, 2 m for the same case as Figure 2. 



7 




Fig. 2 
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